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Abstract 

By considering the one-particle and two-particle scattering data of 
the spin- 1 Heisenberg chain at T = we derive a continuum limit 
relating the spin chain to the ID Bose gas. Applying this limit to the 
quantum transfer matrix approach of the Heisenberg chain we obtain 
expressions for the correlation functions of the Bose gas at arbitrary 
temperatures. 
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1 Introduction 

Bethe's famous investigation on the spectrum of the Heisenberg chain in 1931 [2] 
was performed by a method later termed Bethe's ansatz and was found, after 
suitable modifications, to be applicable to other many-particle systems, e.g. to 
the repulsive Bose gas [18] and the Hubbard model [19]. 

The algebraic reason for the exact solvabihty of these models was found in the 
Yang-Baxter relation [1] and allowed for rather powerful solution techniques on 
the basis of the quantum inverse scattering method [21]. The eigenstates of the 
Heisenberg Hamiltonian are obtained from a trivial vacuum by the action of 
certain operators satisfying the Yang-Baxter algebra, an algebra with quadratic 
relations and structure coefficients given by the i?-matrix of the six-vertex model. 
Within this algebraic Bethe ansatz effective algebraic formulae for scalar prod- 
ucts [22] and norms [16] were derived. 

As a further development, Kitanine et al. [9] were able to express local operators 
of the XXZ chain in terms of the algebra and to calculate static correlation func- 
tions at T = by evaluating scalar products. Thus they rcdcrived the multiple 
integral representations for correlation functions by Jimbo et al. [5,6] and gener- 
alized them to finite longitudinal magnetic fields [8, 10]. 

The next logical step was to introduce finite temperatures with a first appli- 
cation to the generating function [4] of the 2;2;-corrclation. This was achieved 
by Gohmann et al. in [3], where the algebraic structure of [8,22] underlying the 
treatment of ground state correlation functions was combined with the functional 
approach [11-13] for the quantum transfer matrix (QTM). This object describes 
an auxiliary spin chain whose leading eigenvalue is nothing but the partition 
function of the XXZ chain and the eigenvector encodes all thermal correlations. 

In the latest development, with respect to the Heisenberg chain, Sakai [20] suc- 
ceeded in deriving a multiple integral representation for the generating function 
of dynamical correlations for finite temperatures. 

An alternative method to the QTM approach, and actually a more traditional 
treatment of finite temperatures is the so-called thermodynamical Bethe ansatz 
(TBA) introduced by Yang and Yang [23] for the ID Bose gas. Knowing the spec- 
trum in terms of the Bethe ansatz [18], a combinatorial expression for the entropy 
and a minimization of a free energy functional lead to the thermodynamics. For 
the Bose gas this approach yields an integral expression for the thermodynamical 
potential in terms of an auxiliary function to be determined from a non-linear 
integral equation. The same approach, however, for the Heisenberg model yields 
many more integral equations for the bound states of magnons underlying the 
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'string hypothesis'. 

The goal of this paper is the computation of correlation functions for the inte- 
grable Bose gas at finite temperature by use of technical tools developed for the 
spin-1/2 Heisenberg chain. Unfortunately, the direct application of the QTM 
approach fails as it is a concept that is most natural for lattice systems. Here we 
apply these techniques to a lattice system containing the Bose g suitable 
continuum limit. This is where the anisotropic spin-1/2 Heisenberg chain and 
the Bose gas meet. 

In section 2 we review the Bethe ansatz solution of the spin-1/2 Heisenberg chain 
with XXZ anisotropy. In section 3 we identify the continuum limit taking the 
one-particle spectral and the two-particle scattering properties of the Heisenberg 
chain to those of the Bose gas. In section 4 we demonstrate how to rederive 
the famous TBA equations of the Bose gas from thermodynamical equations of 
the Heisenberg chain in the continuum limit. In section 5, a useful generating 
function for correlations is identified and results for this generating function in 
case of the Heisenberg chain are quoted from the literature. These expressions are 
then translated in section 6 into (novel) results for the ID Bose gas by use of the 
continuum limit. This completes our demonstration on how to obtain correlation 
functions for continuum models from those of lattice systems. 



2 The Heisenberg Chain 

The Hamiltonian of the anisotropic XXZ chain on L lattice sites exposed to an 
external longitudinal magnetic field h, 
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is closely related to the transfer matrix of a six- vertex model with local Boltzmann 
weights encoded in the i?-matrix 
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The coupling strength of the Hamiltonian is denoted by J, A = chr/ is the 
anisotropy and u", a — x,y,z denote the usual Pauli matrices. The eigenvalue 
problem of the transfer matrix of the six-vertex model can be solved for instance 
in the framework of the algebraic Bethe ansatz, see e.g. [17]. The energies Em of 
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([T]) and momenta IIm are given by 
d In Am 
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where the transfer matrix eigenvalue Am(A) is obtained from 
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provided that the Bethe equations 
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for the M Bethe roots Xj are satisfied. 



The above hsted Bethe ansatz equations allow for the computation of the entire 
spectrum of the system. For identifying a suitable scaling limit, in which the XXZ 
lattice system turns into the Bose gas in the continuum, it is sufficient to consider 
the one-particle properties and the two-particle scattering data. In particular, we 
are going to identify the ferromagnetic vacuum with the Fock vacuum for Bosons, 
and flipped spins on the ferromagnetic background are considered as Bosons. 
In the one-magnon sector (M = 1) the energy-momentum dispersion simply reads 
(with Eo = -hL/2) 
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and can be considered as some one-particle excitation energy E over the vacuum 
(M = 0) with a chemical potential fi, Ei ~ Eq =: E — fi. 



3 Continuum Limit 



In the desired continuum limit, the ferromagnetic point A = — 1 is approached 
from the region | A | < 1 in a well defined way. Therefore we specify the anisotropy 
of the interaction as A = ch?7 with rj = ivr — ie. Inserting this into the Bethe 
equations as well as the reparametrization A^ = Xj/L of the Bethe roots leads 
to 

'sh(Aj - r]/2)\L _ /I + ieXj/2L^L 
.sh{Xj + 7]/ 2) J ^ \l-iexj/2L 
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where we have used the formula hmn_+oo(l+a;/n)"' = e^ for "LHS" and anticipated 
L ^ oo and exj finite. Introducing a lattice spacing 6 such that the chain has 
the physical length i = 6L, the rescaling exj = iuj reveals the Bethe equations 
to be exactly those of the Bose gas [18], 
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where the repulsion strength can be identified as c = e^/d. Furthermore, the 
single-particle energy-momentum dispersion is 



E = 2J{5vf (10) 

with mass uib = 1/4 (implying J — > oo) and a chemical potential = 2J6^ — h 

directly connected to the magnetic field. For convenience we set = 1/2 as 
Lieb and Liniger [18] already did. 



Parameters for the Spin-| Heisenberg Chain and the Bose Gas 

Let us recapitulate the parameters of the spin chain and the Bose gas with lattice 
spacing 6 and the anisotropy parametrized by A = — cose, e ^ 1: 



XXZ chain 

interaction strength J > 

# lattice sites L, lattice constant 6 

magnetic field h > 

anisotropy A = 5^/2 — 1 



Bose gas 

particle mass itlb = 1/4 
physical length i = L6 
chemical potential fi = 2Je^ — h 
repulsion strength c = 5^/5 



This allows for a well-defined continuum limit for the five lattice parameters while 
keeping the four continuum parameters fixed. In the next section we are going to 
apply the above continuum limit, established at T = 0, to the finite temperature 
formalism of the Heisenberg chain. 



4 Non-linear Integral Equation and Free Energy 

The free energy Fxxz/L per lattice site for the Hamiltonian ([T]) follows within 
the quantum transfer matrix approach [13] as 

'o-^ / . 1 / . , x ln(l + a(a>)) , (11) 
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Figure 1: Depiction of the canonical contour C surrounding the real axis in 
counterclockwise manner within the strip —e/2 < ImA < e/2. 



where the auxiliary function a(A) is calculated from the non-linear integral equa- 
tion 

h 2Jsh^{ie) 



lna(A) = -- + 



T TshAsh(A-ie) 

[ dtu sh(2i£:) , ,^ , ,, 

+ / ^TTx T^TTTx r-ln(l + aH) . (12) 
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Here the temperature is denoted by T and the external magnetic field h in the 
continuum limit (with the fully polarized state in positive direction corresponding 
to the vacuum) takes positive values. In the integrals above, the contour C is a 
rectangular path centered around zero. In the critical regime under consideration, 
with f] = in — ie, its height is restricted by e and the width tends to infinity as 
depicted in figure [H 

Applying the low-T limit (note that J/T — > oo) to the integral equation the 
contribution from the upper line of the contour at = ie/2 vanishes and one can 
restrict both the free argument A and the integration variable u to the lower part 
of the remaining contour. A shift of this line to imaginary part —ri/2 without 
crossing any poles of the inhomogeneity results in the line integral 

. , h 2 J sh^ ri 



T Tsh(A-r//2)sh(A + V2) 



, du sh(2ie)ln(l + a(a7-V2)) 
' 27rish(A-cu + ie)sh(A-cu-i£) ^ ' 



over the real axis. Observing that in the continuum limit the inhomogeneity 
of (fT3|) is nothing but the one-particle excitation energy with some chemical 



7 



potential, cf. ([6]) and (fTOj) . the new auxiliary function reads 

,„„M = -^ + y_^_^_^i„(i + „w) . (14) 

For the free energy we proceed in the same way. First, we note that the Bose gas 
grand canonical potential per unit length (p corresponds to the free energy of the 
Heisenberg chain per lattice constant 6, in detail = {Fxxz / L + h/2) / 5. Second, 
starting from the integral representation ( ITTi) at low temperatures only the lower 
part of the contour contributes. Shifting this integration line to imaginary part 
—7]/ 2 without crossing poles from the kernel we get in the continuum limit 

0=-T / ^ln(l+HH) . (15) 

Looking up the TBA results [23] of Yang and Yang provides the same integral 
representations (fT4|) and (fT5|) describing the thermodynamics. 



5 Finite Temperature Correlation Functions 

The repulsive ID Bose gas is captured by the quantum non-linear Schrodinger 
equation. The Hamiltonian for a system of the physical length I then reads 

Hb= f dz\d,ij\z)d,ij{z)+cij\z)ij\z)ij{z)ij{z)- ^iij\z)ij{z)\ , (16) 

where the operators ip'^ and ip are canonical Bose fields and the mass of the Bosons 
is set to = 1/2 for convenience. Aiming at the number of particles in the 
interval [0,x] the corresponding operator Qi{x) has the explicit representation 

Qi(x) = j dzij\z)i;{z) (17) 



and the ground state expectation value of exp((/?(5i(a;)) can be used to compute 
the density-density correlation function [4] reading 



(^t(a;)^(x)^t(o)^(0)) = i^^(e^QiW) 



Recently Kitanine et al. [7] derived a multiple-integral representation for the 
density-density correlation function for the Bose gas at T = - similar to the 
generating function exp(y9Qi „) of the zz-correlation for the XXZ chain. By iden- 
tifying the operator Qi{x) with Qi^m involving m consecutive sites the definition 
for the spin chain reads 

^ m 
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We already showed [3] that within the quantum transfer matrix approach the 
thermal expectation value (e'^'^i ™)^- has the multiple integral representation 
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and G{ujj, Zk), where G{X, z) is defined by the linear integral equation 
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The integration path F surrounds the origin in counterclockwise manner and has 
to be enclosed by the canonical contour C occurring already in f[T^ . 



6 Density-Density Correlations of the Bose Gas 

For applying the Bose limit to the integral representation (!20l) . let us note that 
the contour C can be chosen to consist of two horizontal lines with imaginary 
parts iie/2, and F similarly, but inside C. 

Proceeding with the low-T limit of (120]) . only the lower integration lines of C 
remain due to the factors a/(l + a). As the upper line of the path F is then no 
longer bounded by the contour C it can be shifted to imaginary part ie/2. Now in 
f l20|) . the integral over this line vanishes in the continuum limit as the factor 



sh(zj + r]) Y" / Vj — ic/2 



sh Zj J \ + ic/2 



x/S 



(23) 
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in the integrand becomes rapidly oscillating asm = x/5^oo for fixed x. 

Like in section HI a simultaneous shift of all remaining (lower) integration lines 
to imaginary part —ri/2 results in 
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for the Bose limit. Suitably reparametrising the variables of integration we find 
the final form 
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As above, we set the physical length of the m consecutive sites to be a; = m6 and 
the limit 6 —>■ led to m ^ oo and to the infinite series on the RHS of the last 
equation. The matrix 
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follows from and as a direct consequence of (12^ the density p{w,p) is the 
solution of the linear integral equation 



2np{u,p) 



{v — p){u — p — ic) 
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(27) 



By use of means rather different from those of [7] we found a new way to derive 
correlations of the Bose gas with even wider applicability, namely to the finite 
temperature case. 



10 



7 Conclusion 

On the level of the Bethe ansatz equations for T = we showed how to relate the 
XXZ chain near the ferromagnetic point to the ID repulsive Bose gas. Applying 
this to the QTM approach reproduces the known TEA thermodynamics of Bose 
particles. Then, and most importantly, we were able to derive a multiple integral 
representation for the generating function of the density-density correlations of 
the Bose gas, valid for arbitrary temperatures. 

Unfortunately, the multiple integral representation is an infinite series with an 
infinite number of integrals to be calculated. In this respect it resembles the Pred- 
holm determinant representation [14,15] but avoids dual quantum fields. Apart 
from the difficulties lying in the infinite series the derivation of explicit results on 
the temperature dependence of the correlation lengths is feasible by considering 
the next-to-leading eigenvalues of the quantum transfer matrix. 

Furthermore it is interesting how to interpret and how to treat other correlation 
functions of the XXZ chain in the considered continuum limit. 
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